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Exchangeable  Phosphorus  Pools  and 
Equilibrium  Characteristics  for  River 
Sediment  as  a  Function  of  Particle  Size 

by  William  F.  James 


PURPOSE:  The  System-Wide  Water  Resources  Program-Nutrient  Sub-Model  (SWWRP-NSM) 
represents  a  library  of  algorithms  for  simulating  nutrient  cycling,  transformation,  and  flux  in 
terrestrial  and  aquatic  systems.  One  feature  of  SWWRP-NSM  is  the  capability  of  simulating 
phosphorus  (P)  equilibrium  fluxes  between  exchangeable  particulate  and  soluble  P  pools  in  the 
aquatic  water  column  as  a  function  of  particle  size  class.  The  objectives  of  this  research  were  to 
quantify  exchangeable  particulate  phosphorus  (P)  pools  and  equilibrium  with  soluble  P  as  a 
function  of  river  sediment  particle  size  distribution  for  use  in  initializing  model  parameters  for 
simulating  P  adsorption  and  desorption  in  river  systems. 

BACKGROUND:  In  agricultural  watersheds,  fertilizer  application  in  excess  of  crop  P  uptake 
requirements  can  result  in  the  enrichment  of  soil  P  levels  over  time  (Sharpley  et  al.  1984;  Zhang 
et  al.  1995,  2004;  Bennett  et  al.  2001).  Numerous  studies  have  shown  that  P  concentration 
increases  with  decreasing  soil  particle  size,  due  to  greater  surface  area  per  volume  and  higher 
composition  of  iron  (Fe)  and  aluminum  (Al)  compounds  that  provide  sorption  sites  for  P  in  finer 
silt  and  clay  soils  (Day  et  al.  1987;  Agbenin  and  Tiessen  1995;  Sinaj  et  al.  1997).  In  contrast,  P 
content  is  much  lower  in  sands  and  associated  primarily  with  calcium  minerals  (Gburek  et  al. 
2005).  Field  runoff  of  particulate  P  during  storms  and  transport  to  receiving  tributaries  is 
governed  by  particle  size  and  density,  as  sandier  soils  tend  to  become  redeposited  in  the  field 
while  the  finest  particles  remain  suspended  in  overland  flow.  The  combined  result  of  these  ero¬ 
sion  and  transport  processes  is  enrichment  of  suspended  particulate  loads  with  P  due  to  move¬ 
ment  of  finer  silts  and  clays  off  the  landscape  (Sharpley  1985). 

During  transport  through  the  river  network,  eroded  particulate  P  can  play  an  important  role  in 
regulating  soluble  P  via  abiotic  processes  such  as  P  equilibrium  between  particulate  and  soluble 
phases  (Froelich  1988).  Differential  sedimentation  due  to  particle  size  and  density  differences 
can  lead  to  transport  of  predominantly  fine-grained  sediment  exhibiting  a  relatively  high  exchan¬ 
geable  particulate  P  concentration.  Equilibrium  processes  can,  thus,  result  in  high  soluble  P  that 
influences  algal  productivity  and  eutrophication,  particularly  as  particulate  loads  are  discharged 
into  impounded  river  reaches  and  coastal  marine  environments.  River  network  models  need  to 
consider  particle  size  distribution  during  transport  in  order  to  better  predict  P  biogeochemical 
cycling  and  influences  on  biota. 

The  SWWRP-NSM  considers  in-stream  processes  for  modeling  P  dynamics,  including  algo¬ 
rithms  describing  P  adsorption  and  desorption  for  a  range  of  particle  size  classes  (Johnson  et  al. 
2008,  2009).  Model  parameterization  requires  partitioning  of  P  into  soluble  and  adsorbed  frac¬ 
tions  as, 
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where,  Psoi  and  Ppart  are  the  soluble  P  and  particulate  exchangeable  P  fractions  (as  a  fractional 
percentage),  respectively,  Cs  is  the  suspended  sediment  concentration  (M  L°),  and  Kd  is  the 

*31 

linear  partition  coefficient  (L  M'  ),  for  particle  size  class  n.  Ppart  and  Psoi  are  directly  related  to 
Kd  according  to  the  equation, 

Ppart  =Kd'Psol  (3) 


For  sands  and  coarse-grained  silts,  we 
would  expect  Ppart  and  Psoi  at  equilibrium 
(i.e.,  no  net  adsorption  or  desorption)  to  be 
relatively  low  due  to  minimal  binding  sites 
for  P  (Figure  1).  Kd  values  should  be  low 
as  well  for  these  particle  size  classes 
because  P  binding  is  weak,  resulting  in 
poor  buffering  between  soluble  and  parti¬ 
culate  phases  due  to  P  disequilibrium.  In 
contrast,  finer-grained  silts  and  clays 
should  exhibit  much  higher  concentrations 
of  Ppart  and  Psoi  at  equilibrium  due  to 
greater  P  binding  potential.  Kd  is  expected 
to  be  much  higher  for  finer  p-grained  par¬ 
ticles,  reflecting  greater  P  buffering  as  a 
result  of  P  disequilibrium  (Figure  1). 
Because  particle  size  distribution  in  river 
suspended  sediment  loads  is  dynamic 
during  transport,  models  are  needed  that 
weight  P  equilibrium  processes  with 
respect  to  particle  size  class,  density,  and 
concentration  in  order  to  more  accurately 
predict  soluble  P.  The  objectives  of  this 
research  were  to  quantify  exchangeable 
particulate  P  pools  and  equilibrium  char¬ 
acteristics  as  a  function  of  river  sediment 
particle  size  distribution  for  use  in  initia¬ 
lizing  model  parameters  for  SWWRP- 
NSM. 


Figure  1.  Conceptual  diagram  showing  relationships 
between  exchangeable  particulate 
phosphorus  (Pparf)  and  soluble  phosphorus 
(Psoi)  for  coarse-  and  fine-grained  particles. 
The  slopes  of  the  red  and  blue  lines 
represent  the  linear  partition  coefficient  (Kd). 
The  dotted  line  depicts  zero  net  P 
adsorption  onto  or  desorption  from  Ppart.  The 
intercept  of  the  lines  with  the  y-axis 
represents  the  concentration  of  the  native 
Ppart  pool  (mg  kg'1).  The  Psoi  concentration 
at  the  point  where  net  P  sorption  is  equal  to 
zero  (i.e.,  intersection  of  the  red  or  blue  line 
with  the  dotted  line)  is  the  equilibrium 
phosphorus  concentration  (EPC).  The  EPC 
is  estimated  as  PpartlKd. 
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METHODS:  Deposited  sediment  was  collected  from  the  Minnesota  River  approximately 
10  miles  upstream  of  its  confluence  with  the  Mississippi  River  at  St.  Paul,  Minnesota,  in  October 
2008.  In  the  laboratory,  the  sediment  was  diluted  (10:1)  with  tap  water  (groundwater)  amended 
with  P  to  a  concentration  of  0.120  mg  L'1  to  maintain  in  situ  P  equilibrium  conditions  (reported 
in  James  and  Larson  2008).  The  solution  was  gently  mixed  on  a  shaker  table  at  20  °C  in  a  dar¬ 
kened  environmental  chamber  before  initiation  of  separation  into  various  particle  size  class 
ranges  (Table  1).  Sieving  and  pipette  techniques  (Plumb  1981)  were  used  to  separate  sediment 
into  particles  >  63  p,  <  63  p,  <  31.2  p,  15.6  p,  <  7.8  p,  and  <  4  p.  Further  separation  was 
achieved  in  a  sedimentation  chamber  by  allowing  particles  to  settle  through  a  100  cm  water  col¬ 
umn  at  20  °C  (Figure  2).  The  sedimentation  chamber  consisted 
of  7  cm  diameter  acrylic  tubing  that  was  attached  to  a  ball  valve. 

A  glass  sedimentation  jar  was  placed  at  the  bottom  of  the  cham¬ 
ber  and  secured  with  a  rubber  stopper  to  seal  the  system.  Par¬ 
ticles  were  gently  shaken  for  24  h  in  P  amended  tap  water  prior 
to  sedimentation.  With  the  valve  closed,  the  sediment  mixture 
was  poured  into  the  upper  portion  of  the  chamber  and  gently 
mixed.  The  valve  was  then  opened  and  sediment  was  allowed  to 
settle  for  various  time  periods,  based  on  fall  velocities  that  were 
computed  using  Stokes  Law.  Finer  particles  remaining  sus¬ 
pended  in  the  chamber  at  the  end  of  the  desired  time  period  were  removed  by  pulling  a  stopper 
located  above  the  collection  jar  and  allowing  the  water  to  drain  into  a  basin.  The  drain  hole  was 
small  (1  cm  diameter)  to  minimize  turbulence  above  the  sedimentation  jar.  Settled  particles  were 
centrifuged,  poisoned  with  chloroform  to  inhibit  biological  activity  (0.1  mL  1%  solution  to 
50  mL  sample),  and  stored  in  a  refrigerator  at  4  °C  until  analysis. 


Table  1.  Particle  Size 
Classes 


>  63  |j 


31.2  to  63  jj 


15.6  to  31 .2  |j 


7.8  to  15.6  |j 


3.9  to  7.8  p 


<  3.9  lj 


Figure  2. 


Sedimentation  chamber  used  for  separating  river 
sediment  into  various  particle  size  classes. 
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Sequential  fractionation  of  particulate  P  was  conducted  according  to  Hjieltjes  and  Lijklema 
(1980),  Psenner  and  Puckso  (1988),  and  Niirnberg  (1988)  to  determine  ammonium-chloride- 
extractable  P  (i.e.,  loosely-bound  P),  bicarbonate  dithionite-extractable  P  (i.e.,  iron-bound  P), 
sodium  hydroxide-extractable  P  (i.e.,  aluminum-bound  P),  and  hydrochloric  acid-extractable  P 
(i.e.,  calcium-bound  P).  A  subsample  of  the  sodium  hydroxide  extract  was  digested  with  potas¬ 
sium  persulfate  to  detennine  nonreactive  sodium  hydroxide-extractable  P  (Psenner  and  Puckso 
1988).  Organic  P  was  calculated  as  the  difference  between  reactive  and  nonreactive  sodium 
hydroxide-extractable  P.  Each  extraction  was  adjusted  to  pH  7  and  analyzed  for  soluble  reactive 
P  (SRP)  using  the  ascorbic  acid  method  (American  Public  Health  Association  2005).  A  portion 
of  loosely-bound,  iron-bound,  and,  to  a  lesser  extent,  the  aluminum-bound  P  fractions  are 
believed  to  be  associated  with  equilibrium  P  exchanges  with  the  water  column.  Total  P  was  cal¬ 
culated  as  the  sum  of  all  measured  extractable  P  forms.  Total  iron  (Fe)  and  aluminum  (Al)  were 
measured  using  atomic  absorption  spectrometry  after  digestion  with  nitric  and  hydrochloric  acid 
(American  Public  Health  Association  2005). 

Aliquots  (-500  mg  L'1  dry  weight  equivalent  to  simulate  suspended  sediment  concentrations  in 
the  lower  Minnesota  River)  of  the  concentrated  suspended  sediment  isolated  for  each  particle 
size  class  were  also  subjected  to  phosphate  (as  KH2PO4)  standards  ranging  from  0  to  1.0  mg  L'1 
for  determination  of  P  adsorption  and  desorption  over  a  24  hour  period.  Tap  water  was  used 
because  it  was  low  in  phosphate  concentration  and  similar  in  ionic  strength  to  Minnesota  River 
water.  Assay  tubes  containing  sediment,  tap  water,  and  known  concentrations  of  phosphate  were 
shaken  unifonnly  in  a  darkened  environment  then  filtered  and  analyzed  for  SRP.  The  sediment 
slurries  were  maintained  under  oxic  conditions  at  a  pH  of  -  8.0  to  8.3  and  a  temperature  of 
-  20  °C  during  the  shaking  and  equilibration  process. 

The  change  in  SRP  mass  (i.e.,  initial  SRP  minus  final  SRP;  mg)  was  divided  by  the  dry  mass 
equivalent  of  sediment  to  determine  the  mass  of  P  desorbed  or  adsorbed  (i.e.,  exchangeable  par¬ 
ticulate  P,  P part',  mg  P  kg'1  dry  mass  sediment).  These  data  were  plotted  as  a  function  of  the  equi¬ 
librium  SRP  (i.e.,  the  final  SRPO  concentration  after  shaking  for  24  hours)  to  detennine  the 
linear  partition  coefficient  (K/,  L  kg'1)  and  the  equilibrium  P  concentration  (EPC)  under  condi¬ 
tions  where  net  sorption  was  zero  (i.e.,  the  Psoi  at  the  EPC;  Froelich  1988).  The  K4  and  Psoi  at  the 
EPC  were  calculated  via  regression  analysis  (Statistical  Analysis  System  1994)  from  linear  rela¬ 
tionships  between  final  SRP  concentration  and  the  quantity  of  P  adsorbed  or  desorbed  at  low 
equilibrium  concentrations.  The  intercept  of  the  regression  equation  represented  Ppart ■  E)ata  were 
also  fitted  to  a  Langmuir  regression  model  using  a  spreadsheet  developed  by  Bolster  and  Horn- 
berger  (2007)  to  estimate  the  Pparl  sorption  maximum  ( Ppart-max )•  The  degree  of  Ppart  saturation 
(DPS)  was  calculated  as  Ppart  divided  by  Ppart-max- 

RESULTS  AND  DISCUSSION:  Concentrations  of  total  P,  loosely-bound  P,  iron-bound  P, 
aluminum-bound,  and  organic  P  increased  with  decreasing  particle  size  class  (Figure  3).  In  con¬ 
trast,  calcium-bound  P  concentrations  were  relatively  constant  over  the  range  of  particle  sizes. 
Total  P  was  composed  primarily  of  calcium-bound  P  for  larger  particle  size  classes  (i.e., 
>  31.2  p,  Figure  4).  This  percentage  declined  with  increasing  particle  size  because  the  calcium- 
bound  P  fraction  remained  constant  while  the  loosely-bound,  iron-bound,  and  aluminum-bound  P 
concentration  increased  as  particle  size  class  decreased.  For  other  P  fractions,  percent 
composition  generally  increased  with  decreasing  particle  size  class.  Iron-bound  P  accounted  for 
the  greatest  percentage  of  total  P  as  particle  size  class  decreased  followed  by  aluminum-bound 


4 


ERDC  TN-SWWRP-1 0-5 
October  2010 

P  >  loosely-bound  P  >  organic  P.  Total  Fe  and  A1  concentrations  also  increased  with  decreasing 
particle  size  class  (Figure  5).  These  trends  suggested  linkages  between  metal  hydroxides  in  river 
sediment  and  extractable  P  forms  involved  in  P  equilibrium  reactions  (Gburek  et  ah,  2005). 


Figure  3.  Variations  in  loosely-bound,  iron-bound,  aluminum-bound,  calcium-bound, 
organic,  and  total  phosphorus  (P)  as  a  function  of  particle  size  class. 
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Figure  4.  Variations  in  the  percentage  of  total  phosphorus  (P)  that  is  loosely-bound,  iron-bound, 
aluminum-bound,  calcium-bound,  and  organic  P  as  a  function  of  particle  size  class. 
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Variations  in  adsorption  and  desorption 
from  the  Ppart  pool  versus  Psoi  for  vari¬ 
ous  particle  size  classes  are  shown  in 
Figure  6.  Generally,  desorption  from 
Ppart  occurred  at  low  Psoi  and  adsorp¬ 
tion  to  PPa,t  occurred  as  Psoi  increased 
above  0.2  mg  L'1.  The  magnitudes  of 
Ppart  desorption  and  adsorption  also 
increased  as  a  function  of  decreasing 
particle  size  class.  Similar  to  iron-  and 
aluminum-bound  P  patterns,  Ppart 
increased  significantly  with  decreasing 
particle  size  class,  indicating  that  the 
exchangeable  P  pool  was  considerably 
larger  for  fine  silts  and  clays  versus 
sands  (Figure  7).  For  particles  less  than 
63  p,  the  Psoi  at  the  EPC  was  high  and 
comparable  to  suspended  sediment  P 
equilibrium  characteristics  in  agricultu¬ 
rally-dominated  systems  (James  and 
Larson  2008).  The  Psoi  at  the  EPC  was 
-0.09  mg  L'1  for  particles  ranging 
between  3 1 .2  p  and  63  p  and  increased 
to  -0.23  mg  L'1  for  particles  less  than 
4  p  (Figure  7).  The  Ppart-max  and  the 
DPS  tended  to  be  inversely  related  to 
particle  size  class  (Figure  8).  These  pat¬ 
terns  may  be  related  to  greater  concen¬ 
trations  of  Fe  and  A1  on  smaller  par¬ 
ticles  with  increased  binding  sites  for  P 
adsorption.  The  Kd  also  increased  as  a 
function  of  decreasing  particle  size, 
suggesting  greater  buffering  of  P 
between  solid  and  aqueous  phases  for 
finer-  versus  coarser-grained  particles 
(Figure  7). 


Figure  5.  Variations  in  total  iron  and  aluminum  as  a 
function  of  particle  size  class. 
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Figure  6.  Relationships  between  mean  (n=3)  exchangeable  particulate  phosphorus  ( Ppart )  and  soluble 
phosphorus  (Pso/)  for  various  particle  size  classes.  Dotted  lines  depict  zero  net  phosphorus 
adsorption  or  desorption.  Vertical  lines  represent  plus  or  minus  1  standard  error  of  the  mean. 
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Figure  7.  Variations  in  mean  (n  =  3)  exchangeable 
particulate  phosphorus  (Pparf),  soluble 
phosphorus  (Psoi),  and  the  linear  partition 
coefficient  (Kd)  as  a  function  of  particle 
size  class.  Vertical  lines  represent  plus  or 
minus  1  standard  error  of  the  mean. 
Different  letters  represent  significant 
differences  (p  <  0.05)  in  the  mean  based 
on  ANOVA  (Duncan-Waller;  Statistical 
Analysis  System  1994). 


Figure  8.  Variations  in  the  theoretical  adsorption 
capacity  for  exchangeable  particulate 
phosphorus  ( Ppart-max )  and  the  degree 
of  phosphorus  saturation  (DPS)  for 
the  exchangeable  particulate 
phosphorus  pool  as  a  function  of 
particle  size  class. 
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Implied  from  these  results  is  the  important  role  that  smaller  particle  size  classes  can  play  in 
regulating  soluble  P  concentrations  via  equilibrium  processes  in  rivers  during  transport.  Clays 
and  finer  silts  have  a  greater  reactive  surface  area  and  higher  Fe  and  A1  concentration  than  sands 
and,  therefore,  greater  potential  for  adsorbing  P.  In  the  watershed,  finer  particles  enriched  with 
adsorbed  P  can  be  preferentially  transported  to  receiving  tributaries  during  runoff  due  to  lower 
density  and  settling  velocity  (Sharpley  1985).  In  large  river  systems,  particle-density  size  differ¬ 
ences  result  in  the  deposition  of  sands  and  coarse  silts  and  transport  of  finer  silts  and  clays  during 
high  flow  events,  which  can  influence  soluble  P  concentrations.  For  the  Minnesota  River,  rela¬ 
tively  high  concentrations  of  Ppart  and  P  fractions  associated  with  Fe  and  A1  for  finer  particles 
were  associated  with  a  Psoi  at  the  EPC  that  exceeded  0.10  mg  L'1.  James  and  Larson  (2008) 
found  that  the  Psoi  at  the  EPC  for  suspended  sediment  in  the  Minnesota  River  was  -0.12  mg  L'1 
during  periods  of  runoff  and  high  flow,  which  is  consistent  with  results  reported  here.  Uptake  of 
soluble  P  loads  generated  by  equilibrium  processes  can  lead  to  eutrophication  and  high  phyto¬ 
plankton  bloom  potential,  particularly  as  loads  enter  impounded  river  reaches.  2-D  and  3D  river 
models  that  simulate  river  channel  P  dynamics  can  be  improved  by  considering  P  equilibrium 
reactions  and  mass  transfer  as  a  function  of  particle  size  during  transport  (Johnson  et  al.  2008, 
2009). 

ADDITIONAL  INFORMATION:  This  technical  note  was  prepared  by  William  F.  James, 
Research  Aquatic  Biologist,  Eau  Galle  Aquatic  Ecology  Laboratory,  Environmental  Laboratory, 
U.S.  Army  Engineer  Research  and  Development  Center.  The  study  was  conducted  as  an  activity 
of  the  Nutrient  Algorithms  work  unit  of  the  System-Wide  Water  Resources  Program  (SWWRP). 
For  information  on  SWWRP,  please  consult  https://swwrp.usace.army.mil/  or  contact  the  Pro¬ 
gram  Manager,  Dr.  Steven  L.  Ashby  at  Steven.L.Ashby@usace.army.mil.  This  technical  note 
should  be  cited  as  follows: 

James,  W.  F.  2010.  Exchangeable  phosphorus  pools  and  equilibrium  characteris¬ 
tics  for  river  sediment  as  a  function  of  particle  size.  SWWRP  Technical  Notes 
Collection,  ERDC  TN-SWWRP-1 0-5.  Vicksburg,  MS:  U.S.  Army  Engineer 
Research  and  Development  Center,  https://swwrp.usace.army.mil/ 
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